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pacs — 

Abstract -Using a new approach, we have obtained a formula for calculating the rotation period 
and radius of planets. In the ordinary gravitomagnetism the gravitational spin (5) orbit (L) cou- 
pling, L ■ S oc L 2 , while our model predicts that L ■ S oc jj L 2 , where M and m are the central and 
orbiting masses, respectively. Hence, planets during their evolution exchange L and S until they 
reach a final stability at which MS oc mL. Rotational properties of our planetary system, binary 
pulsars and exoplanets are in agreement with our predictions. The radius (R) and rotational period 
(D) of tidally locked planet at a distance a from its star, are related by, D 2 oc y-^r R 3 and that 
Roc a. 



introduction. — Kepler's laws best describe the dynamics of our planetary system as 
regards to the orbital motion. However, Newton's law of gravitation provided the theoretical 
framework of these laws. In central potential the orbital angular momentum is conserved. 
In polar coordinates, the gravitational force consists of the ordinary attraction gravitational 
force and a repulsive centripetal force. The Newton's law of gravitation has been successful 
in many respect. However, this law fails to account for very minute gradational effect 
like deflection of light by an intervening star, precession of the perihelion of the planetary 
orbit and the gravitational red-shift of light passing a differential gravitational potential. 
Einstein's general theory of gravitation generalizes Newton's theory of gravitational to give a 
full account for all these observed gravitational phenomena. Einstein treats these phenomena 
as arising from the curvature of space. Hence, Einstein's theory has become now the only 
accepted theory of gravitation. The inclusion of energy and momentum of matter (mass) in 
question leads to the curvature of space, while the inclusion of spin leads to torsion in space. 
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Einstein's theory deals with matter of the former case, while Einstcin-Cartan deals with the 
latter case. Thus, Einstein space if torsion free. In classical electrodynamics the spin of a 
particle is a quantum effect with no classical analogue. However, the spin of a gravitating 
object (e.g., planets) is defined as a rotation of an object relative to its center of mass. This 
is expressed as S = Iuj, where I and to are the moment of inertia and angular velocity of the 
rotating object. The spin is a generally a conserved quantity in physics. Besides the spin, an 
object (to) revolving at a distant r around a central mass (M) with speed v is described by its 
orbital angular momentum. This is defined as L = f x mv. This quantity is also conserved, 
except when an external torque is acting on the object. In quantum mechanics the spin 
and angular momentum of a fundamental particle are quantized. No such quantization is 
deemed to exist in gravitation. To incorporate quantum mechanics in gravitation we invoke 
a Planck-like constant characterizing every gravitational system (6j[T] . This would facilitate 
a bridging to quantum gravity that has not yet been uniquely formulated so far. 

The spin and orbital angular momentum may couple to each other as the case in the 
Earth-Moon system. Therefore, neither the spin nor the orbital angular momentum are 
separately conserved. Their sum is always conserved. A similar coupling occurs in atomic 
system. For instance, because of the spin of the electron an existence of such effect is found 
to be present in hydrogen-like atoms. 

Owing to the existing similarities between gravitation and electromagnetism, some analo- 
gies were drawn which led to gravitomagnctism paradigm. It is believed that an effect 
occurring in electromagnetism will have its counter analogue in gravitomagnctism. 

In this paper we formulate the proper spin-orbit coupling in a gravitational system, 
and then deduce a formula for the spin of a gravitating object. This is done by equating 
the spin-orbit coupling energy to the gravitomagnctic energy. The resulting equation re- 
lates the spin of a gravitating object to its orbital angular momentum. While in standard 
gravitomagnetism, the gravitational spin-orbit coupling, L ■ S cx L 2 , while in our model of 
gravitomagnctism predicts that L ■ S oc jj L 2 . This relation suggests a balance equation, 
toL ~ MS. For this reason any orbiting object must spin to be dynamically stable. So 
planets during their course of evolution exchange L and S, but eventually come to a state 
of stability. The larger the planet the larger its spin. Hence, Jupiter spins faster than other 
planets in the solar system. Equivalently, the spin S oc , where G is the gravitational 
constant and v is the orbital velocity. This formula is found to be consistent when applied 
to our planetary system, exoplanetary system, and binary pulsars. Astronomers have dis- 
covered so far more than 800 new giants planets, but couldn't identify all of their radii and 
spin periods. The present formulation helps identify these latter properties. We consider 
here all possibilities to account for the observational derived data for exoplanetary system 
and their consistency. 
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The gravitational spin-orbit coupling. — The spin - orbit interaction resulting 
from an interaction of the electron spin with the magnetic field arising from electron motion 
in hydrogen-like atom is given by 

Uso= 9s dV z § _ kZe 2 £ g 
4m 2 c 2 r 2 dr 2m 2 c 2 r 3 

where V — r e , Z is the atomic number, k is the Coulomb constant, m is the electron 
mass, c is the speed of light, r is the radial distance of the electron from the nucleus, and 
g s = 2 is the gyromagnctic ratio. 

In gravitomagnctism theory, we have shown that [T] , 

vg' s GM 2 f - 

U so- 9m -^^L.S, (2) 

g ' a is the gravitational gyromagnctic ratio, will correspond to an interaction gravitomagnctic 
energy 



U gm — n ___ Q _a ) (3) 



7T GAIL 2 
3 2mc 2 r 3 

However, Einstein's theory of gravitation employing Schwartzchild metric shows that be- 
cause of space curvature a term of 

Ugr = _ GML ' 2 (4) 
2mc 2 r 3 ' 

appears in the total energy of the gravitating object. 

Thus, eq.(3) and (4) very close to each other. This minute difference between the two 
paradigm should be further explored. Notice that the inclusion of energy momentum in 
Einstein relativity equations leads to the space curvature, whereas the inclusion of spin 
would lead to the space torsion. Einstein's general relativity respects the former but not 
the latter case. While Einstein interprets the precession of planets to curvature of space, 
we ascribe it to the interaction of the spin of planets with gravitomagnetic field of the Sun. 
This interpretation complies with atomic description. Consequently, the gravitomagnctism 
would suggest a curved space due to the presence of electric charge that is tantamount to 
presence of (mass) matter that is known to curve space. 

Assuming the spin-orbit coupling as responsible for precession of perihelion of planetary 
orbits, the spin of a planet of mass m orbiting a star of mass M' can be obtained by equating 
cqs.(2) and (3), i.e., spin-orbit interaction energy equals to gravitomagnetic energy, which 
yields 

_ a m 

S= 7i T7 ' ' 5 

cos M 

where L is the orbital angular momentum of the orbiting planet, a = 9 is the angle 
between L and S directions, e is the orbit eccentricity, and M = (M' + m)o is the total 

4 M ~ M' 
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mass of the system. However, the orbital plane of most planets is inclined to their ecliptic 
with an angle (i), so that in this case eq.(5) becomes 

cost/ M 

Equation (5) suggests that for a system of particle each has a mass of m, and an angular 
momentum of Li, one has a center of mass angular momentum 

„ 772) Li , . 

S ~ 7 

so that the spin represents the center of mass of the angular momenta of the system. 

Owing to the apparent analogy between electromagnctism and gravitomagnetism, one 
has 3 

G -> k , GMm -> kZe 2 , (8) 



so that eq.(l) reads 



where g g is gravitational gyromagnetic ratio analogue. It is not known whether g g = 2 or 
not. Hence, equating eqs.(4) and (8) implies that S — -^-g for g g = 2. This is however not 
correct for the planetary system. Equation (8) agrees with eq.(2) only if g g = 7r^ g' s . 

We remark that several authors have considered the gravitational spin-orbit coupling 
comparing it with the atomic analogue (2] [3]. None of them have derived it from first 
principle, or equivalently didn't show how the gravitational spin magnetic moment is related 
to the spin. This is only done in our recent publication [1]. Mashhoon proposed that the 
analogy between gravity and electromagnetism dictates that charge, q — > —2m. Hence, he 
concluded that the gravitational magnetic moment due to spin is related to spin by /i s = — S 
[4]. In our gravitomagnetic theory, this is however related by the relation \i g = M- S. If we 
write eq.(5) as S = (3 L and the total angular momentum J = L + S, then L ■ S = (3 L 2 . 
Applying eq.(7) in eq.(3) dictates that the curvature term lead to a potential (interaction) 
energy in the atomic system 

Uem= - kZe 2 L 2 
3 2m 2 c 2 r 3 

Therefore, one can write the total potential energy for an electron in hydrogen-like atoms 
in an electrically space as 

E _ kZe 2 | L 2 ir kZe 2 L 2 
r 2mr 2 3 2m 2 c 2 r 3 

Comparing eqs.(l) and (9) reveals that L ■ S = ■§ L 2 for an atomic system. Thus, with 
this understanding the space inside an atom is not flat space (Minkowskian) , but follows 



3 For non-hydrogen-like atom e 2 — > (Ze) e, e — v m, and Ze — > M 
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Schwarzschild pattern, where — > 2 ^J? is the electrical Schwarzschild radius. Hence, 
the metric for a spherically symmetric distribution of nuclear matter can be written as 

2kZe 2 ^ dr 2 
1 



ds 2 = c 2 (l - dt 2 - "" 2kZe2 - r 2 {d6 2 + sine 2 d<p 2 ) . (12) 



One can then interpret the Ruthcrford-a deflection as a consequence of the electrical curva- 
ture inside the atom. This is tantamount to deflection of light by the Sun curvature. The 
electric Schwarzschild radius is equal to twice classical electron radius, R s = 2r c = 2 j°5f . 
Similarly, one would expect a photon to be electrically rcdshiftcd in an electrical potential of 
the nucleus by an amount, z = in hydrogen-like atoms. Hence, in an analogous manner 
when light passe near a central charge it will experience a redshift. This can be written as, 
z = aZ (^f-), where Ac = ^ is Compton wavelength of the electron. Now, when r = Ac, 
then z rn = a Z. This case represents a maximal (quantum) redshift. Thus, owing to the 
one-to-one analogy between gravitation and electromagnctism, one can use Einstein's gen- 
eral relativity to describe electromagnetic phenomena, and Maxwell's equation to describe 
gravitational phenomena. 

The planetary spin and radius. — The origin of spin of planets has not been known 
exactly. One can easily determine the orbital angular momentum of a planet. The spin 
of a planet however requires knowledge of the planet mass, radius, its rotation period and 
its mass distribution inside the planet. Since some planets are solid (rocky) and other are 
gaseous, it is not easily to identify precisely their internal structure. The former ones have 
generally higher rotation rate than the latter. However, orbital periods of planets depend 
on their distance from the Sun and the Sun mass only. We provide here a formula for spin 
or rotational period from its orbital motion only. Or equivalently, we relate the spin to the 
orbital angular momentum for the first time in history. 

Equation (6) can be used to express the planetary spin as 

cos iGm 2 

S = a — . (13) 

v cost' 

This is a very interesting and useful formula that can be used to calculate the spin of a planet 
without resort to its rotation period and its radius. Delauney and Flammarion related the 
spin period of a planet to a host of planetary physical characteristics concluding that there 
is a direct relationship between spin period and mean density |10| . Furthermore, Brosche 
noticed that some planets in a similar size range had spin angular momenta, S, that were 
proportional to the squares of their masses, m [9lll0| 



S cx m 2 . (14) 
Equation (13) agree partially with eq.(12). In rotational dynamics the spin of a rigid body 
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(object) planet is defined by 

S = Iuj, bJ=C %i I = \mR 2 , (15) 

where A is the coefficient of inertia, R the planet's radius, and D is the rotational period of 
the planet. Equation (12) and (14) states that the radius of the planet is 

/ a cos i GmDX 1 ^ 2 . . 

R= tt-t—t;—) ■ (16) 



V 2-kX cos 9 v 



Using eqs. (6) and (14) one can write 



D 1 ( R\ M „ /acosiy'l 



.2 



1/2 



P C 2 \a J to ' ° \ \cos6 ) ' (17) 

where P and a are, respectively, the orbital period and the semi-major axis of the orbiting 
planet. Equation (16) can be written as 

An educated guess can relate a to the ellipticity (flattening/oblateness) of the planet, or to 
the eccentricity of the orbit. Using eq.(14) the rotation spin rate can be written as 



cosi / JGMm 2 a(l - e 2 )\ , s 

» = a W XMR 2 " ' (19) 



Equation (18) can be written as 
cosi L 



h = XMR 2 , L = ^fGMm 2 a{l - e 2 ) . (20) 



cos 9 I\ 

The radius of a planet that is tidally locked to its star, i.e., P = D, is given by (see eg.(16)) 

Rt = Cs lwi a - (21) 

Equation (16) can also be written as, for P = D, 

It is of prime interest to mention that a hypothetical satellite that had a circular orbit radius 
equals to the radius of the planet, R, its orbital period P is given by [10] 



P R = 7777 R3 ■ (23) 



( 47r \ i 
GMj 

Flammarion calculated Pr values for Earth, Jupiter, Saturn, Uranus and Neptune, by ex- 
trapolating Keplers Harmonic Law, as applied to their known satellites [TU]. Only some of 
these period arc in agreement with observation. 
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The radius of a black hole is related to its mass, m, by 

2Gm , , 

R=— s -. 24 

Therefore, the gravitational force for such a planet (spinning black hole) is given by 

_ GmM _ 2 , c 4 D 

r2 -C (1-e )^ — , (25) 

where r = a(l — e 2 ). This clearly shows that a spinning black hole will experience a huge 
gravitational force when orbits any central massive object. 

This force is maximum when the planet is tidally-locked to its star, i.e., P — D. Hence, 
one has 

F^*=C 2 (l-e 2 )^. (26) 

4 

This force is of the order of It is however shown that the maximal force in nature is 
defined by F ma x. = jq [HIE]. It also represents the maximum self-gravitating mass. It is 
thus interesting to see that the gravitational force arising from this case is of the same order 
of this maximal force. For a black hole planet of radius R p tidally-locked with a black hole 
star with radius R s , one has 

R P R S = C 2 a 2 . (27) 

This is an interesting relation connecting the two radii of orbiting black that are tidally- 
locked to their semi major axis. Moreover, it is clear that the existence of such a system 
awaits the future astronomical exploration. 

Binary pulsars. — A pulsar is a rotating neutron star that emits intense beam of 
electromagnetic radiation. A binary pulsar consist of a pulsar orbiting a neutron star or 
white dwarf (companion). The very famous binary pulsar is the PSR B1913+16 discovered 
by Hulse and Taylor [8] The orbital and rotation periods are, respectively, 7.75 hours and 
0.059 sec. 

It is apparent from eq.(l) that for binary pulsars where m ~ M, L ~ S. Hence, eq.(16) 
reveals that 

fl~V? a ' (28) 

The pulsar data can be used to validate eq.(17) and to determine the constant a. We have 
recently place a limit on neutron star higher than those postulated by that of the standard 
one [llj . Accordingly, the radius of neutron star is given by 

*"=(W) WJ =8 °- 9 Ue) (km) ' (M) 
where m e is the mass of the electron and M Q is the solar mass. The radii below in Table 
1 are calculated from eq.(29). We consider here the planctay system, Jupiter satellites, and 
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Saturn satellites. The constant C is calculated using eq.(18) and Tabulated in Tables 1, 2, 
3, 4, 5 and 6 [12] . The average value of C are 0.055, 0.089, 0.077, and 0.068, for binary 
pulsars, Jupiter satellites, Saturn satellites, and planetary system, respectively. Notice that 
for exoplanet and asteroids the constant C takes the average values 7 and 0.5, respectively. 
The higher of C for asteroids may be attributed to the uncertainty associated with the 
observational data related to them. Table 6 can be used to identify cxoplancts that are 
tidally locked by comparing the resulting C for a given system with that of the Moon. Since 
our Moon is tidally locked and well-known, we can consider C = 0.041 for tidally locked 
planets. With this value we complete table 7. Assuming the cxoplanctary system is similar 
to our planetary system, we suggest that C = 0.1, as evident from table 2. With this value 
we calculate the day for some exoplanets as shown in table 8. Notice that Mercury is very 
closed to tidally-lockcd system envisaged in Tables 3 and 7. 
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Name 


to(M ) 


M (M ) 


P (day) 


D(s) 


e 


R (km) 


C 


PSR B1913+16 


1.441 


2.828 


0.323 


0.05903 


0.617131 


97.22 


0.060 


PSR B1534+12 


1.333 


2.678 


0.4207 


0.03790 


0.2736767 


93.48 


0.072 


pQD R9197I IIP 
roll DZLZ 1 -f-llVy 


i 


O 71 o 
Z. ( lo 


U.ooOo 


U.UOUD 


U.Dol 


OA 052 


U.Uoo 


PSR J0737-3039A 


1.34 


2.58 


0.102 


0.0227 


0.0878 


93.73 


0.116 


PSR J1518+4904 


1.75 


2.62 


8.634 


0.04093 


0.2494 


107.02 


0.0414 


PSR B2303+46 


1.16 


2.53 


12.340 


1.06637 


0.658369 


87.21 


0.0076 



Table 1: Binary pulsars primary data 



Name 


m(xl0 24 kg) 


P(day) 


D (hr) 


e 


a (au) 


R(km) 


C 


Mercury 


0.3302 


87.969 


1407.6 


0.2056 


57.91 


2439.7 


0.036390953 


Venus 


4.8685 


224.701 


5832.5 


0.0067 


108.21 


6051.8 


0.009877991 


Earth 


5.9736 


365.256 


23.9345 


0.0167 


149.6 


6378.1 


0.13530133 


Mars 


0.64185 


686.98 


24.6229 


0.0935 


227.92 


3396.2 


0.195060979 


Jupiter 


1898.60 


4332.59 


9.925 


0.04093 


778.57 


71492 


0.087422946 


Saturn 


568.46 


10759.22 


10.656 


0.0489 


1433.53 


60268 


0.111249286 


Uranus 


86.832 


30685.40 


17.24 


0.0565 


2872.46 


25559 


0.079987318 


Neptune 


102.43 


60189 


16.11 


0.0457 


24764 


4495.06 


0.066062714 


Pluto 


0.01305 


89866 


153 


0.244671 


5874 


5874 


0.082682797 



Table 2: The planetary system primary data 



Name 


m(xl0 22 kg) 


a ( km) 


P (day) 


D(day) 


e 


R (km) 


C 


Moon 


7.3477 


384400 27.321582 27.321582 0.0549 


1738.14 


0.041 



Table 3: The Moon (tidally locked to the Earth) primary data 



Name 


m(xl0 2U kg) 


a (10 3 km) 


P(day) 


D(day) 


e 


R (km) 


C 


Io 


893.2 


421.6 


1.769138 


1.769138 


0.004 


1821.6 


0.0116 


Europa 


480 


670.9 


3.551181 


3.551181 


0.0101 


1560.8 


0.0085 


Ganymede 


1481.9 


1070.4 


7.1545535 


7.154553 


0.0015 


2631.2 


0.0051 


Callisto 


1075.9 


1882.7 


16.689018 


16.689018 


0.007 


2410.3 


0.0031 


Elara 


0.008 


11740 


259.6528 


0.5 


0.217 


40 


0.0695 


Himalia 


0.095 


11460 


250.5662 


0.4 


0.162 


85 


0.0483 


Metis 


0.001 


128 


0.294779 


0.294779 


0.0002 


20 


0.3959 



Table 4: Jupiter Satellites primary data 



Name 


m(xl0 2U kg) 


a (10 3 km) 


P(day) 


D(day) 


e 


R(km) 


C 


Miranda 


0.66 


129.39 


1.413479 


1.413479 


0.0027 


235.8 


0.1797 


Ariel 


13.5 


191.02 


2.520379 


2.520379 


0.0034 


578.9 


0.0661 


Umbriel 


11.7 


266.3 


4.144177 


4.144177 


0.005 


584.7 


0.0514 


Titania 


35.2 


435.91 


8.705872 


8.705872 


0.0022 


788.9 


0.0244 


Oberon 


30.1 


583.52 


13.463239 


13.463239 


0.0008 


761.4 


0.0191 



Table 5: Saturn Satellites (tidally locked) primary data 
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Name 


m(xl0 iy kg) 


a (au) 


P (llPO,T\ 


D(dav) 




R {kjri) 


C 


Ceres 


Q7 


7£7 


A R 
4.0 


O H7^ 

y.u < o 


u.u < oy 


A Q7 Q 
4o ( .6 


Q 7E; 


Juno 




z 


RRQ 

z.ooy 


4.O0 


7 01 
I .Zl 


n 0^70 
u.zo < y 


1 on 


r on 
o.yu 


vesta 


oU 


z.ouz 


O.UO 


K Q/l O 


n none 


zoo 


A 71 
4. ( 1 


Eugenia 


n fii 


z. ( zl 


A ACi 
4.4y 


o.oyy 


n nc^i 
U.Uool 


1 1 Q 


1Q If? 


Siwa 


0.15 


2.734 


4.51 


18.5 


0.2157 


51.5 


6.70 


Chiron 


0.4 


13.633 


50 


5.9 


0.3801 


90 


8.48 


Haumea 


41.79 


43.335 


285.4 


3.912 


0.18874 


718 


6.11 


Pallas 


3.18 


2.7707 


3.62 


7.8132 


0.231 


261 


10.03 


Eris 


1.62 


2.385 


3.68 


7.14 


0.231 


199.8 


13.18 



Table 6: Asteroids primary data 



Name 


m(Mj) 


M(M ) 


a (au) 


R(Rj) 


P(day) 


C 


Kepler-34(AB) b 


0.22 


2.0687 


1.0896 


0.76 


288.822 


0.0350 


Kepler-9 c 


0.171 


1 


0.225 


0.823 


38.9086 


0.0395 


KOI-55 c 


0.0021 


0.496 


0.0076 


0.078 


0.34289 


0.0395 


HD 97658 b 


0.02 


0.85 


0.0797 


0.262 


9.4957 


0.0329 


GJ 3470 b 


0.044 


0.541 


0.0348 


0.376 


3.33714 


0.0453 


Keplcr-22b 


0.11 


0.97 


0.85 


0.214 


289.9 


0.0425 


Gl 581 g 


0.01 


0.31 


0.14601 


0.0678 


36.652 


0.0429 



Table 7: Tidally-locked Exoplanets primary data 



Name 


m(Mj) 


M(M ) 


a (au) 


P(day) 


R(Rj) 


D(day) 


WASP-10 b 


3.06 


0.71 


0.0371 


3.09276 


1.08 


0.6566 


XO-5 b 


1.077 


0.88 


0.0487 


4.18775 


1.03 


0.4783 


WASP-16 b 


0.855 


1.022 


0.042 


3.1186 


1.008 


0.1891 


KOI-204 b 


1.02 


1.19 


0.0455 


3.24674 


1.24 


0.1564 


XO-2 b 


0.62 


0.98 


0.0369 


2.61584 


0.973 


0.0993 


TrES-1 


0.761 


0.88 


0.0393 


3.03007 


1.099 


0.1398 


WASP-1 b 


0.86 


1.24 


0.0382 


2.51995 


1.484 


0.0483 


HAT-P-17 b 


0.534 


0.857 


0.0882 


10.3385 


1.01 


2.051 


WASP-55 b 


0.57 


1.01 


0.0533 


4.46563 


1.3 


0.1768 


WASP-6 b 


0.503 


0.888 


0.0421 


3.36101 


1.224 


0.0940 


55 Cnc e 


0.0263 


0.905 


0.0156 


0.736546 


0.194 


0.00578 


OGLE2-TR-L9 b 


4.34 


1.52 


0.0308 


2.48553 


1.614 


0.1079 


OGLE-TR-10 b 


0.68 


1.18 


0.04162 


3.10129 


1.72 


0.04368 


XO-3 b 


11.79 


1.213 


0.0454 


3.19152 


1.217 


1.802 


PSR 1719-14 b 


1 


1.4 


0.0044 


0.0907063 


0.4 


0.000327 


WASP- 14 b 


7.341 


1.211 


0.036 


2.24377 


1.281 


0.4484 


HD 80606 b 


3.94 


0.98 


0.449 


111.436 


0.921 


4444.9 



Table 8: Some non-tidally locked Exoplanets primary data. The day is obtained for C = 0.1 as can 
be guessed from table 2 for planetary system. 
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